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Abstract 

The notion of representative statistical ensembles, correctly representing statistical 
systems, is strictly formulated. This notion allows for a proper description of statis- 
tical systems, avoiding inconsistencies in theory. As an illustration, a Bose-condensed 
system is considered. It is shown that a self-consistent treatment of the latter, using a 
representative ensemble, always yields a conserving and gapless theory. 
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1 Introduction 



Statistical systems are characterized by statistical ensembles. It is crucially important that 
the given statistical system be correctly represented by a statistical ensemble. In other words, 
the chosen statistical ensemble must be representative for the considered statistical system. 
This necessitates a thorough definition of what, actually, a statistical system is, requiring an 
accurate enumeration of all its basic features. The usage of a nonrepresentative ensemble, 
incorrectly representing the considered statistical system, may lead, and often does lead, to 
inconsistencies in the theoretical description of the system. 

The necessity of defining a statistical ensemble that would correctly represent the given 
statistical system was, first, emphasized already by Gibbs [1], who stressed that all additional 
conditions and constraints, imposed on the system, must be taken into account. The prob- 
lem of a proper representation of equilibrium statistical systems by equilibrium statistical 
ensembles was discussed by ter Haar [2,3] and also analized in the review article [4]. 

The aim of the present paper is to formalize the notion of representative statistical en- 
sembles by giving precise mathematical definitions and to generalize this notion for arbitrary 
systems, whether equilibrium or nonequilibrium. The application of the notion is illustrated 
by systems with Bose-Einstein condensate, when the global gauge symmetry is broken. It is 
shown that employing a representative ensemble for a Bose-condensed system results in the 
theory enjoying conservation laws and having no gap in the spectrum of collective excitations. 

Throughout the paper, the system of units will be used with the Planck and Boltzmann 
constants set to unity, h — 1, ks — 1. 

2 Representative Ensembles 

Let us, first, recall several general preliminary definitions that are necessary for precisely 
defining the basic notion of a representative statistical ensemble. 

Physical system is a collection of objects characterized by their typical features distin- 
guishing this collection from other systems. 

For example, a collection of particles can be characterized by their Hamiltonian, that is, 
by their energy operator. 

Statistical system is a many-body physical system, whose typical features are compli- 
mented by all additional constraints and conditions which are necessary for uniquely de- 
scribing the statistical properties of the system. 

Statistical systems are characterized by statistical ensembles. 

Statistical ensemble is a pair {JF, p(t)} composed by the space of microstates T and a 
statistical operator p(t) on that space. 

The space of microstates can be the Fock space or its appropriate subspace. The statis- 
tical operator p(t), generally, is a function of time t. To give p(t) implies to define its form 
p(0) at the initial time t = and to specify the evolution operator U (t) such that 

p(t) = U(t)p(0)U+(t) . (1) 

Therefore, a statistical ensemble can be defined as a triplet 
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The knowledge of a statistical ensemble allows one to find statistical averages. 
Statistical average for an operator A(t) on T is 

< A(t) > = Trrp(t)A(0) = Tr^p(0)i(t) . (2) 

Here the Heisenberg representation of the operator A(t) is assumed, for which 

A(t) = U + (t)A(0)U(t) . (3) 

Representative ensemble is a statistical ensemble equipped with all additional constraints 
and conditions that are necessary for a unique representation of the given statistical system. 

Additional constraints and conditions for statistical systems are usually formulated as 
conditions on statistical averages for some specified condition operators Ci(t), where i = 

1, 2, These operators do not need to be necessarily the integrals of motion, but they are 

supposed to be Hermitian. 

Statistical condition is a prescribed equality for the statistical average of a condition 
operator, 

d(t) = < Ci(t) > = Trp(0)Q(t) . (4) 

Here and in what follows, the trace operation is assumed to be over the appropriate space 
of microstates T . 

Let us consider, first, an equilibrium statistical system, for which the statistical operator 
does not depend on time, 

pit) = p(0) = p. (5) 

The explicit form of the statistical operator follows from the principle of minimal information 
[5]. The latter presumes the conditional maximization of the Gibbs entropy 

,S = -Trplnp (6) 

under the statistical conditions (4), among which one usually distinguishes the definition of 
the internal energy 

E = TrpH (7) 

and the normalization condition 

Trp = 1 . (8) 

The information functional is 

I[p] = S + \ (Trp-l)+f3 (TrpH - E) + (3 £ {jxpC, - C t ) , (9) 

i 

where Ao = InZ — 1 is the Lagrange multiplier preserving the normalization condition (8), 
f3 is the inverse temperature, which is the Lagrange multiplier for condition (7), and flvi are 
the Lagrange multipliers related to statistical conditions (4). 

The minimization of the information functional (9) yields the statistical operator 

P = \ e-? H , (10) 
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corresponding to the grand canonical ensemble with the grand Hamiltonian 

H = H + J2^A. (11) 

i 

The most customary expression for the grand Hamiltonian (11) is 

H = H - fiN , 

where \x is the chemical potential and N is the number-of-particle operator. However, the 
general form of the grand Hamiltonian is given by Eq. (11), in which any condition operators 
can be involved. Thus, an equilibrium representative ensemble is described by the statistical 
operator (10) with the grand Hamiltonian (11). The evolution operator for an equilibrium 
system is 

U(t) = e~ iHt , (12) 
which commutes with the statistical operator (10), because of which 

i j t pit) = [H, pit)] = . 

The general way of obtaining the evolution equations for arbitrary nonequilibrium sys- 
tems is through the extremization of action functionals [6]. In our case, this extremization 
has to be accomplished under the prescribed statistical conditions (4). 

Let the system Hamiltonian be a functional of the field operators ip(x,t) and ip^(x,t), 
that is, H = H[ip], where ip = ip(x,t). The system Lagrangian is 

L[^\ = J if>\x,t) i ^ #M) dx - H[i>\ . (13) 

The action functional, or effective action, under the prescribed statistical conditions (4), 
takes the form 

A[^]= J {LM-5>Ci(f)j dt, (14) 

where Vi are the Lagrange multipliers guaranteeing the validity of the given statistical con- 
ditions. The action functional is defined so that to be a self-adjoint operator, 

A + M = A[ip} . (15) 

Similarly to Eq. (11), the grand Hamiltonian in the Heisenberg representation is 

H[il>]=HW\ + Y,v i C i (t) . (16) 

i 

Then the effective action (14) can be rewritten as 



J y^(x,t)i^(x,t)dx-H[^ dt. (17) 



The extremization of the action functional, requiring that 

8A[1>\ = , 
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with 

yields the evolution equations 



saw _ 0j i^L = . (18) 



6ij>i(x,t) ' <ty(x,t) 

These equations are the Hermitian conjugated forms of each other. 

From Eqs. (17) and (18), it is evident that the evolution equations for the field operators 
can be represented as 

'm^ = W§k (19) 

and its Hermitian conjugated. This should be equivalent to the Heisenberg equation of 
motion 

i-#M) = [#M), H[rl)\\ , 
that is, to the Heisenberg representation for the field operator 

tj}{x,t) = U + (t)^(x,0)U(t) . 
Hence, the evolution operator satisfies the Schrodinger equation 

tj t U(t)==H[^(x,0)]U(t). (20) 

In this way, a nonequilibrium representative ensemble is the set of the given space of 
microstates J 7 , initial statistical operator p(0), and of the evolution operator U(t) defined 
by Eq. (20). An equilibrium representative ensemble is, of course, just a particular case of 
the general nonequilibrium ensemble. 



3 Bose-Condensed Systems 

To illustrate the explicit construction of a representative ensemble, let us consider a system 
with Bose-Einstein condensate. Such systems possess a variety of interesting properties, as 
can be inferred from review works [7-10]. Moreover, theoretical description of these systems 
is known to confront the notorious difficulty of defining a self-consistent approach. The 
theory of Bose-condensed systems is based on the Bogolubov idea [11-14] of breaking the 
global gauge symmetry by means of the famous Bogolubov shift for field operators. The 
condensate wave function, introduced in the course of this shift, has to satisfy the minimum 
of the related thermodynamic potential, which is the stability condition necessary for mak- 
ing the system stable and the theory conserving and self-consistent. At the same time, the 
spectrum of collective excitations, according to the Hugenholtz-Pines theorem [15], has to 
be gapless. The notorious problem is the appearance of the contradiction between the above 
two requirements, when the theory is either nonconserving or gapful. This contradiction 
does not arise only in the lowest orders with respect to particle interactions, when one uses 
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the Bogolubov approximation at low temperatures [11,12] or the quasiclassical approxima- 
tion at high temperatures [16]. However, this contradiction immediately arises as soon as 
the interaction strength is not asymptotically weak and one has to invoke a more elaborate 
approximation. This problem of conserving versus gapless approximations was first empha- 
sized by Hohenberg and Martin [17] and recently covered comprehensively by Andersen [9]. 
The problem is caused by the usage of nonrepresentative ensembles, which renders the sys- 
tem unstable [18]. Here we show that employing a representative ensemble never yields the 
above contradiction, always resulting in a self-consistent theory, being both conserving and 
gapless. 

We consider a system with the Hamiltonian 

+ i J ^(r)^(r')$(r - r')^(r')^(r) drdr' , (21) 

in which the field operators ^(r) = ip(r,t) satisfy the Bose commutation relations, U = 
U(r,t) is an external field, and $(r) = $(— r) is an interaction potential. For describing a 
Bose-condensed system with broken global gauge symmetry, the Bogolubov shift [13,14] has 
to be done through the replacement 

V(r, t) — > V>(r, t) = 77(r, t) + ^i(r, t) , (22) 

where r/(r,t) is the condensate wave function and ipi(r,t) is the field operator of noncon- 
densed particles. The latter field variables are assumed to be orthogonal to each other, 

Jri*(r,t)rl> 1 (T,t)dT = 0. (23) 

It is necessary to emphasize that the Bogolubov shift (22) realizes unitary nonequivalent 
operator representations [18,19]. Accomplishing the Bogolubov shift (22) in Hamiltonian 
(21), as well as in all operators of observables, we get the algebra of observables defined on 
the Fock space generated by the field operators (see details in Refs. [5,18,19]). 

The condensate function is normalized to the number of condensed particles 

N = J \v(r,t)\ 2 dr . (24) 

The Bogolubov shift (22) is only rational when the number of condensed particles (24) is 
macroscopic, which means that the limit 

N 

lim — > 

is not zero, where N is the total number of particles. The latter is given by the average 

N =<N > (25) 

for the number-of-particle operator 

N = [ ^(rMr) dr , (26) 
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in which the Bogolubov shift (22) is again assumed. The statistical averaging in Eq. (25) 
and everywhere below is over the Fock space 

Substituting the Bogolubov shift (22) into Hamiltonian (21) gives in the latter the terms 
linear in ipi, because of which the average < ipi > can be nonzero. This, however, would 
result in the nonconservation of quantum numbers, e.g., of spin or momentum. Therefore, 
one has to impose the constraint for the conservation of quantum numbers, 

<^(r,t) > = 0. (27) 

Defining the self-adjoint condition operator 

A(t) = J [A(r, t)4(r, t) + A*(r, t)fa(r, *)] dr , (28) 

in which A(r, t) is a complex function, we may represent constraint (27) as the quantum 
conservation condition 

< A(t) > = . (29) 

In this way, there are three statistical conditions. The first condition is the normalization 
(24) for the number of condensed particles. Condition (24) can be represented in the standard 
form (4) by defining the operator 

N = ij \v(r,t)\ 2 dr , (30) 

in which 1 is the unity operator in the Fock space J 7 {fa). Then Eq. (24) reduces to the 
statistical condition 

N = < N > . (31) 

The second condition is the normalization (25) for the total number of particles. Equivalently, 
instead of normalization (25), we may consider the normalization condition 

N 1 = <N 1 >, Ni ee J 4{r)fa{r) dr (32) 

for the number of uncondensed particles Aq = N — N . And the third condition is the 
conservation condition (29). Respectively, the effective action, which is now a functional of 
the two field variables, r/{r,t) and fa(r, i), with taking account of the statistical conditions 
(29), (31), and (32), becomes 

fa] = / {L + VoN + HiNi + A) dt. (33) 

Here L = L[ip] is the Lagrangian (13) under the Bogolubov shift (22) and A = A(t) from 
Eq. (28). The quantities fj, , /xi, and A(r, t) are the Lagrange multipliers guaranteeing the 
validity of the corresponding statistical conditions. Introducing the grand Hamiltonian 

H[ V , ^} = H- fi N - HiNi - A , (34) 

in which H = H[ip], with shift (22), and the effective Lagrangian 



Lfo iM = J 



d d 
rf (r, t)i — rj(r,t) + 4 (r, t) i — fa (r, t) 



dr - H[ V , Vi] , (35) 
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for the action functional (33), we get 

A[ V , Vi] = J L[riM dt 



(36) 



The evolution equations follow from the extremization of the action functional (36), that 
is, from the variations 

<L4[77,Vi]_ n (37) 



5r]*(r : t) 







and 



SA[ v ,ij 



H = 0. 



(38) 



These equations, as is clear from Eqs. (35) and (36), are equivalent to the equations of 
motion 

i 1 nfr t) - 5 J^hM m 
% dt V{T,t) ~ 5r]*(v,t) m 



and 



(40) 



dt riV ' ' #I(r,t) 

One has to substitute here Hamiltonian (34) under the Bogolubov shift (22). Accomplishing 
the variation in Eq. (39), we get 



*^(M) 



2m 



+ U-Ho ?7(r) + 



+ j $(r-r') [\ V (r')\ 2 r)(r)+X(r,r')] dr' , (41) 

where the time dependence in the right-hand side, for short, is not explicitly shown, U = 
U(r,t), and the notation for the correlation operator 



X(r,r') = 4(r')Mr'Mr) +^(0^(0^(0+ 
+rf (r>i(r>i(r) + ^WiWiW 
is used. The variation in Eq. (40) gives 



(42) 



2m 



+ U-in) ^i(r) + 



+ /$(r-r') ^(r 7 )! (r) + v%r')v(r)M r ')+v(r')v(r)i>l(r')+X(r,r') dr' . (43) 



To get an equation for the condensate wave function, we have to take the statistical 
average of Eq. (41). For this purpose, we introduce the normal density matrix 



Pl (r,r') = < ^(r>i(r) > , 



(44) 



the anomalous density matrix 



(7i(r,r 



(45) 
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and their diagonal elements, giving the density of noncondensed particles 



pi(r) = px(r, r) = < Vi(r)^i(r) > (46) 

and the anomalous average 

o-i(r) = <ri(r,r) = < Vi(r)^i(r) > . (47) 

The quantity |o"i(r)| can be interpreted as the density of paired particles [19]. The total 
density of particles 

p(r) =p (r)+pi(r) (48) 

consists of the condensate density 

p (r) = |77(r)| 2 (49) 

and the density of noncondensed particles (46). Averaging Eq. (41), we find the equation 
for the condensate wave function 

d / V 2 \ 

i ai'^ t)= [-2^ + u ->"')^ )+ 

+ J $(r - r) p(r')rj(r) + pi(r, r)^(r') + <7i(r, r')rj*(r')+ < ^i(r')^i(r')^i(r) > dr' . (50) 

Equations (43) and (50) are the basic equations of motion for the field variables i](r,t) and 
ipi(r, t). These equations, according to Eqs. (39) and (40), are generated by the variation 
of the grand Hamiltonian (34). The latter, in agreement with Eq. (20), also defines the 
evolution operator U(t), which satisfies the Schrodinger equation 

tj t U(t) = H[ V (r,0), ^(r,0)]f)(t). 
Thus, the representative ensemble for a Bose-condensed system is the triplet 

It is important to stress that the so defined representative ensemble possesses a principal 
feature making it different from the standardly used ensemble having the sole Lagrange 
multiplier fi = fi 1 . But then the normalization condition (24) cannot be guaranteed. Then 
the evolution equation for the condensate wave function is not a result of a variational 
procedure. For an equilibrium system, this means that the number of condensed particles 
No does not provide the minimum of a thermodynamic potential, which implies the system 
instability. All notorious inconsistencies in theory, manifesting themselves in the lack od 
conservation laws or in the appearance of an unphysical gap in the spectrum, are caused by 
the usage of nonrepresentative ensembles. 
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4 Green Functions 

The equations of motion (43) and (50) allow us to derive the evolutional equations for the 
Green functions. To this end, we shall use the compact notation denoting the set {rj,tj} by 
the sole letter j, so that the dependence of functions on the spatial and temporal variables 
looks like 

/(12 ...n) = /(ri, ti, r 2 , t 2 , ■ ■ ■ , r n , t n ) . 
The product of the differentials dvjdtj will be denoted as d(j), so that 

n 

d(12...n) = Y[drj dtj . 

3=1 

We shall employ the Dirac delta function 

5(12) = <J(n - r 2 ) 5(t ± -t 2 ). 

For the interaction potential, we shall use the retarded form 

$(12) = $(n - r 2 )<J(fi - f 2 + 0) . (51) 

The matrix Green function (7(12) = [G aj g(12)] is a 2 x 2 matrix, with a, (3 = 1,2, and 
with the following elements: 

G n (12) = -i< Ty>i(l)Vi(2) > , Gi 2 (12) = -i < 7^(1)^(2) > , 

G 21 (12) = -i< f 4(1)4(2) > , G 22 (12) = -i< f 4(1)^(2) > , (52) 

where T is the time-ordering operator. 
Let us introduce the operator 

the condensate effective potential 

V(12) = 6(12) J $(13)1^(3)1^(3) + $(12)77(1)^(2) , (54) 
and let us rewrite the correlation operator (42) in the form 

1(12) = ^(2)^(2)77(1) + ^(2)77(2)^(1) + 77*(2)^(2)^(1) + ^1(2)^(2)^(1) . (55) 
We also define the matrix correlation function X(123) = [X a p(123)] with the elements: 
X n (123) = - < TX(12)4i(3) > , X 12 (123) = - < f X(12)Vi(3) > , 

X 21 (123) = - < TX + (12)4i(3) > , X 22 (123) = - < f X+(12)^i(3) > . (56) 
From the equations of motion (43) and (50), we find the equations 

G n (12) - J V(13)G ll (32) d(3)- 
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-J $(13) [77(l)77(3)G 2 i(32) +iX n (132)] d(3) = 5(12) , 

(i J- - tf^ G 12 (12) - | \/(13)G 12 (32) d(3)- 

- J $(13) [77(1)77(3)^22(32) + iX 12 (132)] d(3) = , 

(-z ^ - C 21 (12) - I V*(13)G 21 (32) d(3)- 

- J $(13) [77*(l)77*(3)Gn(32) +iX 21 (132)] d(3) = , 

(-i ^ - K^j C 22 (12) - I ^*(13)G 22 (32) d(3)- 

- I $(13) [r/*(l)r/*(3)G 12 (32) + 7X 22 (132)] d(3) = 5(12) . (57) 
The self-energy £(12) = [£ a /j(12)] is a matrix whose elements are defined by the relations 

[Eii(13)Gn(32) + £ 12 (13)G 21 (32)] d(3) = 

= I V(13)G n (32) d(3) + I $(13) [77(1)77(3)G 21 (32) +iX n (132)] d(3) , 

I [£ 11 (13)G 12 (32) + £ 12 (13)G 22 (32)] d(3) = 
= I \/(13)G 12 (32) d(3) + I $(13) [77(1)77(3)^22(32) +iX 12 (132)] d(3) , 

I [£ 21 (13)G 11 (32) + £ 22 (13)G' 21 (32)] d(3) = 
= I V*(13)G 21 (32) d(3) + I $(13) [77*(l)77*(3)G n (32) +iX 21 (132)] d(3) , 

I [£ 21 (13)G 12 (32) + £ 22 (13)G 22 (32)] d(3) = 

= I \/*(13)G 22 (32) d(3) + J $(13) [77*(l)77*(3)Gi 2 (32) +^X 22 (132)] d(3) . (58) 

Let us introduce the matrix condensate propagator C '(12) = [C Q( g(12)], with the elements 

C n (12) = -i77(l)77*(2) , C 12 (12) = -277(1)77(2) , 

C 21 (12) = -277*(1)77*(2) , C 22 (12) = -^(1)77(2) , (59) 
The latter have the properties 

C n (21) = C 22 (12) , Cn(ll) = C 22 (ll) , C 12 (21) = C 12 (12) , 

C 21 (21) = C 21 (12) , C* u (12) = -C 22 (12) , C* 2 (12) = -C 21 (12) . (60) 
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The binary Green function is a matrix 5(123) = [B a/3 (123)], 



5(123) = C(12)G(23) - ip (2)G(13) + X(123) , 



(61) 



whose elements are 



B n (123) = Cn(12)G n (23) + C 11 (22)G' 11 (13) + C 12 (12)G 2 i(23) + X n (123) , 

B 12 (123) = C , i 1 (12)G 12 (23) + Cn(22)Gi 2 (13) + Ci 2 (12)G 22 (23) + X 12 (123) , 

5 21 (123) = C 22 (12)G 2 i(23) + C 22 (22)G 21 (13) + C 21 (12)G 11 (23) + X 21 (123) , 

5 22 (123) = C 22 (12)G' 22 (23) + C 22 (22)G 22 (13) + C 2 i(12)Gi 2 (23) + X 22 (123) , (62) 

and where po(l) = |^(1)| 2 - 

With Eqs. (59) and (61), relations (58), defining the self-energy, can be rewritten in the 
matrix form 



in which the delta function 5(12) in the right-hand side is assumed to be factored with the 




(63) 



Then the equations of motion (57) acquire the matrix representation 




(64) 



unity matrix 1 = [6 a p\ and 



73 



I 







-1 



is a Pauli matrix. 

Introducing the inverse propagator 




(65) 



allows us to transform Eq. (64) into 




(66) 



An equivalent representation, following from Eq. (66), is 




(67) 



For the self-energy, using Eq. (63), we have 




(68) 
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The equations for the Green functions are to be complimented by the equation for the 
condensate wave function (50), which, introducing one more anomalous average 



£(12) = < ^(2)^(2)^(1) > , 



(69) 



can be represented as 



V? 
2m 



+ 



J $(12) [p(2)rj(l) + pi(12)ij(2) + < 7 1 (12V(2) + £(12)] d(2) . (70) 

It is the equations for the Green functions and the equation for the condensate func- 
tion, which become mutually incompatible in the standard approach, while employing the 
representative ensemble renders the theory self-consistent in any approximation. 



5 Theory Self-Consistency 

One usually confronts inconsistencies in theory considering a uniform equilibrium Bose- 
condensed system. Then, in any given approximation, one gets either a nonconserving 
theory, that is, an unstable system, or one finds an unphysical gap in the spectrum, which, 
actually, again corresponds to an unstable system [9,18]. To analyze this problem, we pass 
now to the case of an equilibrium uniform system, when U — 0. 
Then we use the Fourier transform for the Green function 



G(12) = J G(k,u)e i{k - 



(2tt)' ' 
in which 

r 12 = r l — r 2 j ti2 = ti — t 2 ■ 

By their definition in Eq. (52), the Green function elements possess the properties 

G n (21) = G 22 (12) , G 12 (21) = G 12 (12) , G 21 (21) = G 21 (12) . (71) 

Therefore the corresponding Fourier transforms satisfy the relations 

Gn(-k, -u) = G 22 (k, u) , Gi 2 (-k, -co) = G 12 (k, u) , 

G 21 (-k,-u) = G 21 (k,u). (72) 
Assuming that the system is isotropic, one has 

G a/3 (-k,uj) = G a p(k,Lj) (73) 

for all a, j3. Combining Eqs. (72) and (73), we find 

Gn(k, -u) = G 22 (k, uj) , Gi 2 (k, -u) = Gi 2 (k, u) , 
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G 21 (k,-cu) = G 21 (k,cu) . (74) 
Also, for a uniform equilibrium system, one has [14] the equality 

G 21 (k,uj) = G 12 (k,uj) . (75) 

Fourier-transforming the self-energy 

\dkdu 



£(12) = J £(k,cj)e i(k ' ri2 -^ 12) ' 



(27T) 4 

and, similarly, the inverse propagator (65), we have for the latter 

k 2 

G-\k, u) = r 3 cu - — + fi - £(k, u) . (76) 
2m 

Then Eq. (67) reduces to 

G- 1 (k,cj)G(k,cj) = i , (77) 

where i = [8 a p\. 

From Eqs. (76) and (77), it follows that G _1 (k, u), hense, also £(k, u), have the same 
symmetry properties as G(k,u>). In particular, 

£ Q/3 (-k,u;) = Y, a p(k,uj) , £i 2 (k, -u) = T, 12 (k,uj) , £ 2 i(k, -w) = £ 2 i(k, w) , 

E 2 i(k,w) = £ 12 (k,cu) , En(k,-w) = £ 22 (k,cu) . (78) 
The matrix equation (77), explicitly, is the system of equations 

( W ~~ 2~m + /il ~~ Sn ) Gn ~~ Sl2G<21 = 1 ' 
(v - ^ + /ii - £n^j G12 - £i 2 G 22 = , 
^-w - ^ + /ii - £ 22 ^J G 2 i - E 2 iGn = , 

(-oj - ^ + A*i - £ 22 ) G 22 - £ 21 G 12 = 1 , (79) 

where, for short, = G a p(k, uo) and £ Q/3 = £ Q/3 (k, u;). The solutions to these equations 
are 

uo + k 2 /2m + £n(k, a;) — fi± 



Gn(k,w) 



^ /, v £i 2 (k,cu) 

G 12 (k,.) = --^— r , (80) 



with the denominator 



D(k,u) 



2 



k 

u - £n(k,w) +fi! 

2m 



x 
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X 



k 2 

lo + h E 22 (k, w) - Hi 

2m 



+ E? 2 (k,cu). (81) 



The solutions for E 2 i(k, u;) and G 22 (k, a;) are defined by the symmetry properties (72) to 
(75). 

The excitation spectrum is given by the poles of the Green functions, that is, by the zero 
of denominator (81), 

D(k,e k ) = 0. (82) 

Equation (82) can be represented as 

e k = X - [E n (k, e k ) - E 22 (k, e k )\ ± y/d* - Ef 2 (k, e k ) , (83) 

with the notation 

k 2 1 

u k = — + - [E n (k, e k ) + E 22 (k, e k )} - m . (84) 
2m 2 

Denominator (81) enjoys the property 

D(k,-u) = D(k,u) . 

Consequently, if e k is a solution of Eq. (82), then — e k is also its solution, which is in 
agreement with the form of Eq. (83). 

For an equilibrium uniform system, the Bogolubov shift (22) is equivalent to the sepa- 
ration of the zero-momentum term in the expansion of the field operator over plane waves. 
The shift itself has meaning only under the normalization condition (24), in which N ~ N, 
that is, the zero-momentum state is macroscopically occupied. The latter becomes possible 
when the single particle spectrum touches zero. Therefore, the necessary condition for the 
existence of Bose-Einstein condensate is 

lime fc = 0. (85) 

This is to be complimented by the stability condition 

Re e k > , Im e k < . (86) 

This condition should be kept in mind when choosing the sign plus in front of the square 
root in spectrum (83). 

Taking limit (85) for spectrum (83), we notice that, according to properties (78), 

E n (k,0) = E 22 (k,0). (87) 

By using perturbation theory for a stable system, one can show [15] that in all orders of the 
theory 

S Q/3 (0,0)>0. (88) 
Then the necessary condition (85) yields the expression for the chemical potential 

/i 1 = E 11 (0,0)-E 12 (0,0), (89) 
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which is the Hugenholtz-Pines relation [15]. 

On the other hand, we have Eq. (70) for the condensate wave function. For an equilibrium 
uniform system, with no external potential U, all densities do not depend on the spatial and 
temporal variables, 

p (r)=p , pi(r) = pi, <7i(r) = <7i, p(r) = p . (90) 

The condensate wave function reduces to the constant 

V( T , t) = V = VPo ■ (91) 
Then we substitute into Eq. (70) the Fourier transforms for the interaction potential 

for the normal density matrix (44), 

Pi(ri,r 2 )=/n fc e*-^L, 
and for the anomalous density matrix (45), 



Similarly, the Fourier transform for the anomalous average (69) is 



6(ri,r 2 ) = jte*'™ — 



(27T) 

As a result, Eq. (70) gives 



Po 



P% + f(n k + a k + ) $ k 7^ ■ (92) 

Generally, expressions (92) and (89) do not coincide with each other, their difference 
being 

r ( & \ x dk 

Po ~ Pi = P®o + / \n k + a k + —— — — - - 

-E n (0,0) + E 12 (0,0). (93) 

This is the general expression for the difference between the Lagrange multipliers p an d Pi 
for an arbitrary equilibrium uniform Bose-condensed system. 

Usually, one does not distinguish between the Lagrange multipliers p and pi, which 
implies setting p — Pi — > 0. However, as is evident from Eq. (93), there is no any reason 
for requiring that this quantity be zero. As an illustration, we may resort to the Hartree- 
Fock-Bogolubov approximation, in which £ k = and 

n k $ k , E 12 (0, 0) = p $o + J °k®k ■ 
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Relation (89) then yields 



Hi = p% + / {n k - a k )<& k 



dk 



(94) 



(2tt) 3 " 



The difference of the chemical potentials (93) becomes 



Ho - Hi 




(27T) 



3 ' 



(95) 



which, certainly, is nonzero [20]. 

In this way, the introduction of the additional Lagrange multiplier makes the theory 
completely self- consistent. All inconsistencies that often arise in other works, such as the 
appearance of a gap in the spectrum, system instability or a distortion of the phase transition 
order, are caused by neglecting the difference between the multiplier fii an d the multiplier 
Ho- It is worth emphasizing that the introduction of the Lagrange multiplier fj, for preserving 
the normalization condition (24), from the mathematical point of view is strictly necessary. 
In other case, the employed ensemble would not be representative, hence, could not correctly 
describe the Bose-condensed system with broken gauge symmetry. 
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